Influence of the anion potential on the charge ordering in quasi-one dimensional 

charge transfer salts 
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We examine the various instabilities of quarter-filled strongly correlated electronic chains in the 
presence of a coupling to the underlying lattice. To mimic the physics of the (TMTTF^X Bechgaard- 
Fabre salts we also include electrostatic effects of intercalated anions. We show that small displace- 
ments of the anion can stabilize new mixed Charged Density Wave-Bond Order Wave phases in 
which central symmetry centers are suppressed. This finding is discussed in the context of recent 
experiments. We suggest that the recently observed charge ordering is due to a cooperative effect 
between the Coulomb interaction and the coupling of the electronic stacks to the anions. On the 
other hand, the Spin-Peierls instability at lower temperature requires a Peierls-like lattice coupling. 

PACS: 75.10.-b, 75.50.Ee, 71.27. +a, 75.40.Mg 



Site-centered Charge Density Wave (CDW) or Bond 
Order Waves (BOW) states are generic in one- 
dimensional (ID) chains at commensurate fillings. It is 
known that a strong enough short range electronic re- 
pulsion leads to a (4k f) charge instability |l],|j Molecu- 
lar crystals such as quasi-lD charge transfer salts [QJ|] 
present a very rich physics due to the interplay between 
electron-electron and electron-phonon interactions. Sev- 
eral systems have been observed to show transitions to- 
wards charge ordered phases where the molecules of the 
conducting stack exhibit unequal electron densities. The 
sulfur series of the Bechgaard-Fabre salts family, namely 
(TMTTF) 2 X (X = PF 6 , AsF 6 ), consists of quarter-filled 
(n = 1/2) molecular chains with very anisotropic trans- 
port properties E[ . Below the Mott localization crossover 
temperature evidences for a true transition towards a 
symmetry broken state have been recently provided by 
dielectric response measurements Q and NMR || . This 
lower symmetry phase is characterized by the disappear- 
ance of all central symmetry centers in the crystal [0. 
The Spin-Peierls (SP) transition observed at lower tem- 
perature Q is characterized by a freezing of remaining 
spin fluctuations (opening of a spin gap) and it is accom- 
panied by a tetramerization of the chains Q . 

The interplay between electronic correlations and lat- 
tice effects in ID quarter filled chains have been addressed 
in a number of previous theoretical studies ||-|| where in- 
teresting modulated phases have been proposed. Gener- 
ically, lattice modulations (or BOW) were shown to be 
always accompanied by CDW's of weaker amplitudes || . 
However, the potential created by the intercalated an- 
ions have also been invoked as an extrinsic cause to 
the 4hp dimerization. In that weaker lattice cou- 

pling is sufficient to produce a similar 2kp instability 
(tetramerization). In this work, we extend our previous 
model of Ref. to include the crucial effect of the anion 
potential. The following approach is based on the real- 



istic two-dimensional (2D) structure of the charge trans- 
fer salts depicted in Figure (for more details see also 
Figs. 1 & 8 of Ref. §1 or Fig. 3.3 of Ref. §). We only 
focus here on the most anisotropic compounds of the sul- 
fur series where inter-chain charge transfer can be safely 
neglected. 

Our model consists of a 2D array of electronic (ex- 
tended) Hubbard chains (hereafter called E-chains) cou- 
pled to classical lattice degrees of freedom; 

+ V ^ II, L ,»,.., + -ffclas + ^anions i (1) 



where i (j) are site (chain) labels, riij ;a — c[j. a Ci,j; a 
and riij — fiij-^ + We have included a nearest 

neighbor (NN) interaction V. Small local displacements 
of the molecules along the E-chains can couple strongly 
to the electrons via modulations of the single particle 
hoppings of the Peierls type, 



t(2p-l,j)=4(l + 5^_ ld ), 
t(2p,j) = t° 2 (l + gS^), 

with an elastic energy 

H clas = lK B J2(S^) 2 - 



(2) 



(3) 



The electron-lattice couplings have been (partially) ab- 
sorbed in the re-definition of the bond modulations SP 4 so 
that the strength of the lattice coupling scales like 1/ Kb- 
Let us now consider the role of the charged X~ anions lo- 
cated on the so-called A-chains as shown in Fig. |[ They 
introduce a new periodicity of 2 lattice spacings on the 
E-chains which, a priori, might lead to a dimerization of 
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the bare electronic transfer integrals along the chain i.e. 
t? ± t%. Hereafter, the mean value t = /2 is set to 

1 and we assume t ? — t§ = O.lt. However, the qualitative 
conclusions of this work do not really depend whether t\ 
and t® are assumed to be equal or not. For simplicity, 
we shall take g = 1. The anions can play a dominant 
role if they are allowed to undergo small displacements 
(along some arbitrary direction) leading to local changes 
of the on-site electronic energies. In lowest order in these 
displacements 8^ one gets, 

^anions = Y.^ X6 P-J ~ 6 P,j+l) n 2p,j ( 4 ) 
P,3 

where p labels the positions of the anions in the chain 
direction. Note that all displacements (for molecules and 
anions) are defined with respect to the high-symmetry 
(high-T) phase shown in Fig. [p. The parameter A (0 < 
A < 1) accounts for the non-equivalent location of each 
given anion with respect to the two nearest molecules in 
each of the two neighboring chains. Hereafter, results are 
obtained for an intermediate value A = 0.3. 
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FIG. 1. Schematic structure of the (a,b-c) plane (see Fig. 8 
of Ref. used in this paper: electronic chains (grey sites) 
are separated by anionic chains (black sites) with one an- 
ion for two organic molecules. Due to the complete charge 
transfer of one electron per anion the electronic chains are 
quarter-filled on average. The arbitrary displacements (exag- 
gerated for clarity) of the anions with respect to symmetric 
positions are shown by arrows. Central symmetry centers C, 
C, Ci and C2 in the high-symmetry configuration are indi- 
cated on the plot. 

Following the method of Ref. || based on Exact 
Diagonalisation (ED) of small clusters supplemented 
by a self-consistent procedure we obtain the lowest 
(T=0) energy lattice configuration without any assump- 
tion on the super-cell order of the broken symme- 
try ground state (GS). The condition that the energy 
E({Sfj,Spj,}) is minimum with respect to the sets of 
distortions/displacements {Sfj} and {Spji} reads, 

K B 8f d + t(4- ld . lt ,Cij., a + h.c.) = , 



KaS^j, + \((n 2p ,j>) - (n 2 p+ij'-i>) (5) 
+ (ri2p+i,j>) - (nzpj'-i) = 0. 

Here (...) is the GS mean value obtained by ED (using the 
Lanczos algorithm) of Hamiltonian ([!]) on periodic M xL 
clusters. Typically, M = 2 units ( 2 electronic chains and 
2 anionic chains) were considered in the transverse direc- 
tion (so that modulated structures with k± = or 7r 
can spontaneously appear) with up to L = 12 (L/2 = 6) 
electronic (anionic) sites in the E-chains (A-chains). The 
non-linear coupled equation (|B|) can be solved by a gen- 
eralization of the iterative procedure of Ref. jlO) . Note 
that since spin and charge are conserved in each elec- 
tronic chain, each of them can be diagonalized indepen- 
dently. Note also that, for finite lattice coupling, finite 
size effects are in general quite small. 

Our numerical results show a very rich phase diagram 
with mixed CDW-BOW structures. Quite generally the 
largest supercell unit is four lattice spacings long along 
the E-chains and has two basic (E-chain+A-chain) units 
in the transverse direction. Therefore, in order to discuss 
our results, it is convenient to parameterize the real space 
behavior of the various observables in the following way; 

7T 

Ariij = p 4kF (-iy cos (k±j) + p2k F cos(-i + fc ± j + $ 2 fc F ) , 

Sfj = ^ kF (-lT cos (k^j) + Sg F cos (li +k± j+^ kF ) , 

Sfa - 5$ cos (k±j) + 5*{-\y cos (fcjj) , (6) 

where Anij = {(ni,j) —n)/n is the relative change of the 
electronic density. Note that the periodicity of the 4kp 
(resp. 2kp) modulation corresponds to 2 (resp. 4) lattice 
spacing of the E-chains while a uniform (hj_ = 0) or a 
staggered (fcj_ = 7r) arrangement occur in the transverse 
direction. Note that both the 4fcp-BOW (dimerization), 
the 2fc F -BOW sequence X-O-X-0 ($f fcp = 0) and the 
2fc_F-CDW sequence B-B-B-B ($2fc F = tt/4) preserve at 
least one of the central symmetries of Fig. [y (X = —A). 

Before discussing the generic phase diagram of the 
model, it is interesting to comment on the special 
limit Ka = 00, i.e corresponding to fixed anion posi- 
tions (Sp j =0). In that case, the E-chains decouple 
from each other and the problem reduces to the single 
Peierls-Hubbard chain The phase diagram of this 

model H exhibits, for sufficiently large lattice coupling, 
an instability towards two different translation symme- 
try breaking phases D\ and D2] (i) in the weak cou- 
pling regime (let's say U/t < 3), a strong 2k F BOW with 
&2k F = 7r / 4 i- e - corresponding to a X-X-X-X type of 
sequence of the bonds occurs. This modulation coexists 
with a weaker 2kp site-centered CDW (A-O-A-0 type of 
sequence of the on-site relative charge densities) , and an 
even weaker Akp (i.e. A-A-A-A) CDW component; (ii) 
at larger U/t, a superposition of a lattice dimerization 
(4kp BOW) together with a tetramerization (2kp BOW 
with <&2k F = 0) both coexisting with a weak 2/cf CDW 
component (with $2fe F = tt/4) occurs. It is interesting to 
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notice that the dimerization observed in experiments on 
Bechgaard-Fabre salts Q| is absent in the weak coupling 
D\ phase. In addition, the D2 phase, a natural can- 
didate to describe the low temperature SP phase, still 
exhibits one center of symmetry (e.g. C2) while recent 
experiments [Q point towards the disappearance of all 
center of symmetries below the charge ordering transi- 
tion. Therefore, we believe that the anions may play an 
important role in stabilizing new lower symmetry CDW- 
BOW phases. 
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FIG. 2. Phase diagram versus the two electron-lattice cou- 
plings for U — 6, V = 2 and A = 0.3. Dashed regions are 
beyond the validity of the model. The different phases are 
depicted in Fig. 0. U denotes the uniform phase (with only 
bare dimerization). The numerical values of the order param- 
eters along the dashed lines indicated by the arrows can be 
found in Fig. 

The T=0 phase diagram of our model is shown in Fig. ^ 
for typical parameters in Bechgaard-Fabre salts. A rich 
variety of charge ordered phases are found and schemat- 
ically depicted in Figs. [|. We note that (at least on 
finite clusters) finite values of the electron-lattice cou- 
plings are necessary to stabilize these translation sym- 
metry breaking phases. At small Peierls coupling, the 
coupling to the anion displacement field generates a 
k e = (4&; F ,0) electronic CDW simultaneously with a 
dimerization in the chain direction (hereafter named as 
"a-phase"). This charge ordering is accompanied by a 
uniform = (0, 0) displacement of the anions with re- 
spect to their symmetric positions as shown in Fig. ||(a). 
On the other hand, the D2 phase of the single chain 
system (see Ref. Q) at large enough 1/Kb coupling is 
weakly affected by the coupling to the anions and con- 
tinuously evolves into a two-dimensional in-phase array 
(i.e. kx = 0) of mixed 2fc F -CDW/2fc F -BOW/4£; F -BOW 
chains (named as b-phase here) which extends in a wide 
region of the phase diagram. Note that, as soon as 1/K a 
is finite, this mixed CDW-BOW phase is accompanied 



by a k^4 = (7T.0) component of the anion displacement 
field (see Fig. 0(b)) so that the translation symmetry is 
fully preserved in the transverse direction. As for the 
single chain D2 phase, the two-dimensional b-phase ar- 
rangement preserves the symmetry centers C2 and C . 
Novel interesting phases are obtained in the region of in- 
termediate 1/K a and 1/Kb couplings. The c-phase of 
Fig. ||(c) [resp. d-phase of Fig. ||(d)] is realized by a 
superposition of a uniform k = (0, 0) displacement and 
a k^ = (71-, 0) [resp. k^ = (7r,7r)] modulation of the 
anion positions. The charge ordering is characterized 
by the superposition of a 4fc F -CDW with a 2fc F -CDW 
whose phase < $2/c F < 7r/4 depends on the model 
parameters. Similarly the lattice modulation contains a 
4fc F -component (dimerization) superposed with an out- 
of-phase (i.e. < <S> 2 3 kp < tt/4) 2fc F -BOW. While all the 
E-chains are in-phase in the c-phase, they are alternat- 
ing in the d-phase. It is important to notice that only 
the uniform phase and the b-phase preserve at least one 
center of symmetry of the crystal. 
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FIG. 3. Schematic representation of the various phases 
found in this work. Different average charge densities are in- 
dicated by different levels of grey, (a) and (b) [resp. (c) and 
(d)] have two [resp. four] non-equivalent sites. While (a) and 
(b) show two and three non-equivalent bonds respectively, (c) 
and (d) have four. 

The modulation amplitudes are plotted vs 1/Kb in 
Figs. ||(a,b) for two characteristic values of the electron- 
anion coupling constants. These plots show clearly that 
the Peierls coupling stabilizes BOW's while the coupling 
to the anions rather stabilizes 4fc F -CDW. We believe that 
the competition between these two effects produces the 
interesting mixed 4-component CDW-BOW phases at in- 
termediate 1/K a and 1/Kb values. In the c-phase and 
d-phase the phase of the 2fc F CDW [resp. 2fcp BOW] in- 
creases rapidly from ~ to ~ 7r/4 [resp. decreases from 
~ 7r/4 to ~ 0] for increasing Peierls coupling to eventually 
become locked to the value 7r/4 [resp. 0] in the b-phase 
(analogous to the ID D2 phase). Note that all center 
of symmetries are removed by the 4fc F -CDW component 
in the a-phase, c-phase and d-phase and/or also by the 
2kp CDW and/or BOW components in the c-phase and 



3 



d-phase as soon as $2/c F ^ 7r/4 and/or $f fe ^ 0. Note 
also that electronic interaction (e.g. U) tends to suppress 
both 2kp CDW and BOW components as seen in Ref. |§ . 



0.8 



2k F BOW 
_ 4k F BOW 
at 06 O-O 2k F CDW 
*> 4k F CDW 



1/K A =1 



(a) 



£ 0.4 
Q. 

i a-phase 

0.2 



0.0 



<0 06 
<D 

~o 

3 

£ 0.4 
Q. 

E 

" 0.2 



□ 



& 0-= 



b-phase 



0.0 



1/K A =1.47 








(b) 










■ 


a-phase 




/ 












M"' 


b-phase 


a- 


--n- 








# • 


fl O 








e- - e — 



0.4 



0.8 

1/K„ 



1.2 



1.6 



FIG. 4. Amplitudes of the various 2&F and 4&f CDW and 
BOW components (as indicated on the plot) as a function of 
the Peierls coupling 1/Kb, for U = 6, V = 2, A = 0.3 and 
two different anion-electron coupling constants as indicated 
on the plots. The correspondance to the ordered phases of 
Fig. ^ is shown. 

We finish by a discussion of our results in the light 
of the recent experiments |^,|| providing evidences of a 
charge modulated state in (TMTTF) 2 X (X=PF 6 , AsF 6 ). 
Although the spontaneously dimerized phase S (L>2 
phase) of the (uniform) single chain Peierls-Hubbard 
model is a reasonable candidate for the charge ordered 
state [Q, it seems in contradiction with the conclusion 
of the dielectric response measurements due to its re- 
maining central symmetry center. Hence, although the 
single chain Peierls-Hubbard model contains most of the 
relevant features of the experimental system, a correct 
understanding still requires to consider the effective cou- 
pling of these (electronically decoupled) chains via the 
anion potential. We believe that the dimerized a-phase 
stabilized above a critical value of 1/ K a might be a fair 
description of the ordered phase with 2 inequivalent sites 
in agreement with NMR spectroscopy results || . Clearly 
Coulomb interaction (V term) also cooperates with the 
electrostatic anion potential to generate the 4/c^ CDW 
instability. The small additional 2k F BOW and CDW 
components of the c-phase (or d-phase) might also ex- 
plain the additional translation symmetry breaking (and 
the spin gap opening) occurring at the SP transition. 
Indeed, although at T = we obtain first order lines be- 



tween the various GS (by varying the parameters), con- 
tinuous transitions could appear at finite temperature 
(and once lattice dynamics is introduced) between let say 
the uniform phase and the a-phase at high temperature 
and the a-phase and the c-phase (or d-phase) at lower 
temperature. Incidently, this scenario then implies that 
the charge ordering must be accompanied by a (probably 
quite small) global shift of the anions with respect to their 
symmetric locations and that the subsequent SP transi- 
tion is linked to an additional (even smaller) = (tt, 0) 
or k.A = (tt, tt) modulation. 
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